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Can lightning be a noise source 
for a spherical gravitational wave antenna? 
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The detection of gravitational waves is a very active research field at the moment. In Brazil the 
gravitational wave detector is called Mario SCHENBERG. Due to its high sensitivity it is necessary 
to model mathematically all known noise sources so that digital filters can be developed that maxi¬ 
mize the signal-to-noise ratio. One of the noise sources that must be considered are the disturbances 
caused by electromagnetic pulses due to lightning close to the experiment. Such disturbances may 
influence the vibrations of the antenna’s normal modes and mask possible gravitational wave signals. 
In this work we model the interaction between lightning and SCHENBERG antenna and calculate 
the intensity of the noise due to a close lightning stroke in the detected signal. We find that the 
noise generated does not disturb the experiment significantly. 
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I. SCHENBERG’S FEATURES 

The detection of gravitational waves is presently an ac¬ 
tive research field and there are several detectors either 
being built or in operation around the world, involving 
countries like USA, Japan, Germany, The Netherlands, 
Italy and Australia. In Brazil the gravitational wave de¬ 
tector is called Mario SCHENBERG. This resonant-mass 
detector will operate in coincidence with some of the de¬ 
tectors in other countries. It is now in an advanced pro¬ 
cess of installation at the University of Sao Paulo, in the 
capital of Sao Paulo State. It is expected that SCHEN¬ 
BERG will be tested with three of its nine transducers 
in 2006. 

Resonant-mass gravitational wave detectors like 
SCHENBERG detect the presence of a signal basically 
through the vibration of its first quadrupolar modes. 
Transducers strategically positioned on the antenna’s 
surface are able to transform these vibrations into elec¬ 
tric signals that can be digitally processed. As one can 
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see from Figure ^ SCHENBERG’s antenna is spherical; 
it is made of a CuAl alloy (6% Al) and has 65 cm in di¬ 
ameter. More details on its functioning can be found in 
reference [Q. 

According to General Relativity the quadrupolar 
modes are expected to be excited by gravitational waves 
p]. Due to its characteristics in SCHENBERG the first 
of these modes are tuned to « 3.2 kHz. However, if an 
electromagnetic pulse generated by a lightning stroke hits 
the antenna most of its normal modes should be excited, 
including the quadrupolar ones. Therefore lightning is 
expected to create noise in the detector perhaps masking 
detectable gravitational waves. 

In a previous work we have modelled the interaction 
between cosmic rays and SCHENBERG |jj. In that pa¬ 
per the interaction was basically thermodynamical. In 
this work the approach is different: the electromagnetic 
wave generated by the lightning stroke hits the antenna 
and transfers momentum to it. As a first approximation 
to the calculations we will assume that SCHENBERG 
consists only of the spherical antenna, ignoring for now 
the cryogenic chambers. 

In what follows, after a brief description of how 
SCHENBERG would respond to the influence of an 
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electromagnetic wave we will present the mathematical 
model for lightning strokes. Then the lightning force den¬ 
sity is calculated since it will be used in the model of the 
interaction between the detector and the lightning. Re¬ 
sults and conclusions are the closing sections. 


In the case of SCHENBERG 6 ~ 1.5 mm. Thus an ap¬ 
preciable fraction of absorption loss should happen prac¬ 
tically at the antenna’s surface. 


B. Reflection loss in the near field 


II. ELECTROMAGNETIC WAVES IN 
CONDUCTING MEDIA 

A lightning creates an electromagnetic pulse that prop¬ 
agates with the speed of light. We want to investigate 
what happens when this pulse reaches the antenna, which 
is a conducting solid medium. For the description of 
the behavior of electromagnetic waves in conducting me¬ 
dia one can use results found in the literature 0. Two 
types of losses are encountered by an electromagnetic 
wave striking a metallic surface: 1) the wave is partially 
reflected from the surface (reflection loss), and 2) the 
transmitted portion is attenuated as it passes through 
the medium. This latter effect is called absorption loss. 
Absorption loss is the same in either the near or the far 
field and for electric or magnetic fields. Reflection loss, 
on the other hand, is dependent on the type of field and 
the wave impedance. 

The transition region between the near and the far 
fields is around A/27T, where A is the source’s wavelength. 
In the case of SCHENBERG, the relevant wavelength is 
A = c/v o = 3 x 10 8 /3175 ~ 9.4 x 10 4 m, thus implying 
a transition region around 15 km. We will consider dis¬ 
tances smaller than 15 km because they include stronger 
strokes. Within this distance the far field limit would be 
suited only for higher frequencies, which are much weaker 
then the lower ones, as one can see from Figure [3 For 
these reasons we conclude that the electromagnetic fields 
generated by lightning that hit the ground close to the 
detector are in the near field limit. 


A. Absorption loss 


It can be shown that when an electromagnetic wave 
reaches a good conductor it is reduced to 1/e = 0.369 
of its initial amplitude after travelling the distance (all 
units in this paper are in the MKS system) 



The distance S is called skin depth. It varies with the 
frequency of the wave, as can be seen in Table [I] 0. The 
symbols p and er refer, respectively to the permeability 
(equal to 47r x 10 - ' H/m for free space) and to the con¬ 
ductivity (equal to 5.82 x 10 7 mhos/m for copper). 

Then for a medium with thickness t the absorption loss 
is given by 


dB. 


In the near field the electric (E) and the magnetic ( H ) 
fields must be considered separately, since the ratio of the 
two is not constant. In this limit reflection loss varies 
with the wave impedance, defined by Z = E/H. For 
the same source frequency a high-impedance (electric) 
field has a higher reflection loss than a low-impedance 
(magnetic) field. 

In the case of SCHENBERG, mostly made of copper 
and sensitive to frequencies around 3.2 kHz, we can ex¬ 
pect a minimum reflection loss around 50 dB 0. In fact, 
for low frequency plane waves, as the ones considered 
here, reflection loss accounts for most of the attenuation. 


III. MODEL FOR THE LIGHTNING 

The phenomenon of lightning is still not fully under¬ 
stood but several models have already been proposed to 
explain it [!]. In brief, a cloud-to-ground lightning flash 
is usually composed of several intermittent discharges 
called strokes, which are made up of a leader phase and 
a return stroke phase. The leader initiates the return 
stroke, which is an upward travelling wave. A relatively 
large electric field exists in the leader-return stroke chan¬ 
nel and this field produces ionization and results in a 
current. The power input renders the channel very lumi¬ 
nous and causes its rapid expansion, producing thunder. 
More details on the phenomenon can be found in 7]. 

In our analysis we will consider the effect of a first 
return stroke on SCHENBERG. The idealized stroke is 
drawn in Figure 0 adapted from 

The electric field intensity E and the magnetic flux 
intensity B produced by the channel current i(z,t ) at 
a distance D from the vertical channel of height H are 
given by (gj 
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V. THE MODEL OF THE INTERACTION 
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where eo is the permittivity, po the permeability of free 
space and c is the speed of light in vacuum. 

There are several models for the channel current i(z, t) 
all trying to describe actual strokes as accurately as 
possible. Since we are interested mostly in the intensity 
of the disturbance caused by a stroke in the detection of 
gravitational wave process we will use a simple expression 
for the current, which accounts for the basic features of 
a typical lightning. 

We assume that the stroke has a typical length of 5 
km. It takes 10/xs to reach the maximum current of 30 
kA and lasts typically 50/rs. With these values we model 
the current as a function of time as a Poisson distribution: 


The mechanical model of a spherical gravitational wave 
antenna has already been presented in the literature. We 
will follow here the notation on m for the case that there 
are no transducers coupled to the antenna and that only 
source of force on the sphere is the e-m wave. Then the 
Fourier transform of the radial amplitude of motion of 
the detector surface under the influence of this wave is 
given by 


dm(^) — k I m 



i ' W ° i 2 

+ »q' w + w 0 


/ M, (5) 


where /(w) is the Fourier transform of the force density, 
the constant k is defined by 

k = a{R)R 2 / (gM eff N) 


and 


Im, — 


7r 7r 

J j Y 2 m (#, <t>) sin 2 0 si n <t> d9 d<fi. 


9=0 d>=0 


i(t) = ate xt , (3) 

with a = 8.155 x 10 9 A/s and A = 10 5 s -1 . See Figure 
Q] for the plot of this function. We will adopt here the 
Bruce-Golde model 0, for which i{z,t) = i(0,t). Also, 
we will assume that the lightning stroke is quite close to 
the detector, at a distance D = 1500 m, which ensures 
a high amplitude to the electromagnetic front that will 
disturb the antenna. 


The list of constants used in the equations is presented in 
Table HU They are based on recent data from SCHEN- 
BERG [ll|. The functions Y 2m (#, <j>) are the spherical 
harmonics and can be found in textbooks [1] . As usual, 
wo = 27ri'o. The limits of integration are such that they 
cover the half-sphere that is in the direction of the light¬ 
ning stroke. 


VI. RESULTS 


IV. DETERMINATION OF THE LIGHTNING’S 
FORCE DENSITY 

From the expressions CD. © and © we were able 
to compute the Poynting vector for the electromagnetic 
(e-m) radiation that reaches SCHENBERG using the for¬ 
mula S (t) = /x ( )' 1 E(t) xB(f). The coordinate system used 
is depicted in Figure |3J with D fixed at 1.5 km. From 
this expression we found the density of momentum of 
the wave, p (t) = S (t)/c 2 . The plots of the electric and 
magnetic fields at this close distance are shown in Figure 
0 

By differentiating p(t) we obtained the expression for 
force density of the e-m wave, f(f), which for the distance 
we chose has modulus 

f(t) = 2.21 x 10“ 3 e" 10 ’ 5 * (t - 10” 5 ) + 


+e 2 xio 5 t (_ 2 .83 x io 3 t 2 + 0.768 t - 2.42 x 10~ 7 ) (4) 
with units of VTs 2 m~ 3 in the MKSA system. 


From equation © we could determine the amplitude 
of the motion caused by the e-m wave ( a m (t )) in the five 
degenerate quadrupole modes of SCHENBERG, given by 
the index m = —2, —1,0,1, 2. We found that modes m=l 
and m=-l are not disturbed at all by this wave. Modes 
m=2 and m=-2 are equally disturbed, slightly more than 
mode m=0. 

In gravitational wave detection the adimensional am¬ 
plitude h m = a m /R is preferred as a measure of the an¬ 
tenna’s motion. In terms of this figure the results we 
obtained on the motion of the modes with time is shown 
in Figure [U] Notice that the modes present a maximum 
adimensional amplitude around 3.4 x 10 -24 . 

VII. CONCLUSIONS 

The gravitational wave detector SCHENBERG is ex¬ 
pected to reach a maximum sensitivity of h ~ 10 -2 ° for 
impulsive waves when fully operational in a few years. 
The noise caused by a typical, close lightning stroke on 
the detector was found in this work to be approximately 
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three orders of magnitude smaller than this expected sen¬ 
sitivity. 

Had the electromagnetic shielding due to the metallic 
cryogenic chambers be taken into account in the calcu¬ 
lations this noise would be even less significant. There¬ 
fore we conclude that it is unlikely that lightning strokes 
should cause detectable noise while SCHENBERG is run¬ 
ning. 

This result is relieving since it may happen that 
SCHENBERG does not run in coincidence with other 
detectors from time to time, and it must be free of as 
much kinds of noise as possible. The region in which this 
detector is been built is prone to summer thunderstorms 
during approximately 2 months, accompanied by many 
lightning strokes. We have just shown that such strokes 
should not directly disturb the data significantly. Occa¬ 
sionally electric fluctuations may happen in power supply 
due to lighting strokes but these can be easily ruled out 


with the continuous monitoring of the experiment. 

Although the maximum adimensional amplitude due 
to lighting strokes found above is larger than the ampli¬ 
tudes due to many monochromatic gravitational waves 
pj this kind of noise should not disturb the detection of 
this kind of waves since the signal is integrated in time. 
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LIST OF FIGURE CAPTIONS: 

Figure ^ Schematics of the SCHENBERG detec¬ 
tor. The spherical antenna is carefully suspended inside 
metallic chambers that keep it in vacuum at 4 K. 

Figure [2 Plot of the Fourier transform of the light¬ 
ning’s current. 

Figured Schematics of an idealized return stroke with 
the geometrical and physical parameters used in the cal¬ 
culations. 

Figure 0J Plot of the behavior of the current of the 
lightning as a function of time, modelled according to a 
Poisson distribution. 

Figure [3J Plots of the electric (a) and magnetic (b) 
fields as functions of time for D = 1500 m. 

Figure m Adimensional envelop of amplitudes as a 
function of time for modes (a) m=0 and (b) m=2. Notice 
that the amplitude decays slowly due to the antenna’s 
high mechanical quality factor. 
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FIG. 1: 




|l(f)| (A/hz) 
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Frequency 

5 

5 

(Hz) 

for Cu for A1 

60 

8.51 

10.90 

100 

6.60 

8.46 

1 k 

2.08 

2.67 

3.2 k 

1.63 

1.49 

10 k 

0.66 

0.84 

100 k 

0.20 

0.28 

1 M 

0.08 

0.08 

10 M 

0.02 

0.02 


TABLE I: Skin depth S for copper and aluminum, in mm. 


Name 

Symbol 

Value 

Radius at 4K 

R 

0.3239 m 

Geometric constant 

a(R) 

2.862 

Density at 4K 

Q 

8077.5 kg/m 3 

Effective mass 

M e ff 

288 kg 

Normalization constant 

N 

0.142kg 3 

Quadrupole mode frequency 

vo 

3175 Hz 

Mechanical quality factor 

Q 

2 x 10 7 


TABLE II: Values of the constants used in the calculations. 


Ligthning Current in frequency domain 



FIG. 2: 









